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We develop an exact method for computing Casimir forces and the power of radiative heat transfer 
between two arbitrary nanostructured surfaces out of thermal equilibrium. The method is based 
on a generalization of the scattering approach recently used in investigations on the Casimir effect. 
Analogously to the equilibrium case, we find that also out of thermal equilibrium the shape and 
composition of the surfaces enter only through their scattering matrices. The expressions derived 
provide exact results in terms of the scattering matrices of the intervening surfaces. 
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I. INTRODUCTION 

In recent years impressive technical advances have 
prompted intensive experimental and theoretical investi- 
gations of proximity phenomena, originating from quan- 
tum and thermal fluctuations of the electromagnetic 
field existing in the vicinity of all bodies. In gen- 
eral, these phenomena can be grouped in two broad 
classes, namely equilibrium phenomena one side, and 
non-equilibrium phenomena on the other. Well-known 
examples of equilibrium phenomena are the Casimir ef- 
fect, and the Casimir-Polder atom- wall forces On 
the other hand, a much studied non-equilibrium phe- 
nomenon is provided by non-contact radiative heat trans- 
fer between two closely spaced bodies Q. In recent 
times much interest has been devoted to the new field of 
Casimir and Casimir-Polder forces between macroscopic 
bodies and/or atoms out of thermal equilibrium. These 
phenomena are intensely investigated now both theoreti- 
cally 0, [1[ and experimentally [y| . We should also like to 
mention the very interesting phenomenon of non-contact 
quantum friction [7J. The problems of heat transfer, 
Casimir forces and quantum friction in a system of two 
plane-parallel plates at different temperatures, in relative 
uniform motion in a direction parallel to the plates, have 
also been investigated recently ||. 

While the material dependence of the above phenom- 
ena has been well studied in simple planar geometries, 
there exists presently much interest in exploring more 
complicated geometries. Indeed the highly non-trivial 
geometry dependence of the near-field opens up the pos- 
sibility of tailoring the features of the radiation field for 
new applications, that range from micro- and nanoma- 
chines operated by the Casimir force [Toll . to designing 
the thermal emission of photonic crystals [ill] . For nearly 
flat surfaces, the shape dependence can be studied us- 
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ing the so-called proximity force approximation (PFA), 
which amounts to averaging the plane-parallel result over 
the slowly varying distance between the opposing sur- 
faces of the bodies. It is well known however that the 
PFA has only a limited range of validity, and it can lead 
to inaccurate predictions for surfaces with deep corruga- 
tions [12J. For this reason, it is widely recognized today 
that more accurate methods are needed to describe ar- 
bitrary geometries. For systems in equilibrium, much 
progress has been made recently. New powerful numer- 
ical techniques to compute the Casimir force, based on 
the Green's function approach and on the path-integral 
approach have been reported [H, 03] . Another approach 
that is being vigorously developed is based on the mul- 
tiple scattering formalism, that was first introduced long 
ago [l5| to study the Casimir energy for a system of per- 
fect conductors of arbitrary shape. The scattering ap- 
proach is actually closely related to the Green's function 
method, and we address the reader to Ref. for fur- 
ther details on the connection between the two methods. 
There exist today several variants of the scattering ap- 
proach, that have been developed to deal also with real 
materials (for a comparative review see [ill]): one of the 
variants [It], is better suited for dealing with compact 
objects not too close to each other, and it is based on a 
multipolar expansion of the e.m. field. Another variant 
[l9| is instead better adapted to deal with planar-like 
structures in close proximity, and it uses a decomposi- 
tion of the e.m. field into plane waves. Contrasted with 
these important advances in the mathematical techniques 
for computing equilibrium Casimir forces, the theory is 
much less developed for problems out of thermal equilib- 
rium. In fact, the Casimir force out of thermal equilib- 
rium has been investigated only in the plane-parallel case 
0, while the problem of computing the near-field radia- 
tive transfer between two spheres was addressed only very 
recently [201 ] . using a dyadic Green's function approach 
(see also the recent experiment |21[). Investigating in 
depth the shape-dependence of proximity effects out of 
thermal equilibrium is indeed very interesting, in view of 
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potential applications, because out of equilibrium there 
exists a richness of behaviors, associated for example with 
resonances in the spectrum of surface excitations, that 
are absent at equilibrium [H, 0| • 

In this paper we develop an exact method for com- 
puting Casimir forces and the power of heat transfer be- 
tween two arbitrary plates out of thermal equilibrium. 
The method is based on a generalization of the scatter- 
ing approach, that has proven so successful in equilib- 
rium Casimir problems. We consider the variant of the 
scattering approach [lj| that is best suited for planar-like 
nanostructured surfaces at close separations, like those of 
11, 12]. In this paper, we shall only present the deriva- 
tion of the basic formulae, leaving concrete numerical 
applications for a successive exposition. The main new 
result is the demonstration that also out of thermal equi- 
librium, the shape and material dependence enter only 
through the scattering matrices of the bodies involved, 
analogously to what has been found for systems in ther- 
mal equilibrium. We remark that our results provide 
exact expressions in terms of the scattering matrices of 
the intervening bodies. Of course, the scattering ma- 
trix is in principle a complicated object, but there exist 
methods, both analytical and numerical, for computing 
it accurately. Many new geometries have been considered 
recently in Casimir investigations, and the corresponding 
scattering matrices have been estimated for this purpose. 
The formulae derived in this paper permit to consider 
these geometries out of thermal equilibrium. 

The plan of the paper is as follows: in Sec. 2 we present 
the general derivation of the correlators for the e.m. field 
in the gap between two arbitrarily shaped plates at dif- 
ferent temperatures. In Sec. 3 the correlators derived 
in Sec. 2 are used to obtain an exact expression for the 
Casimir force and the power of radiative heat transfer 
between the plates. Finally, in Sec. 4 we present our 
conclusions and outline directions for future work. 



II. GENERAL PRINCIPLES: RYTOV'S 
THEORY. 

We consider a geometry of the type considered in Refs. 
i.e. a cavity consisting of two large plates plates 
at temperatures T\ and T 2 , the whole system being in 
a stationary configuration. The shapes of their oppos- 
ing surfaces can be arbitrary, apart from the assumption 
(usually implicit in scatte ring approaches to equilibrium 
Casimir problems [ll], EH Il9j l) that there must exist be- 
tween them a vacuum gap of thickness a > bounded 
by two parallel planes. This condition excludes from our 
consideration interpenetrating surfaces, like the one stud- 
ied in the second of Refs. [14]]. We also assume for sim- 
plicity that the plates are thick, in such a way that no 
radiation from outside can enter the gap. 

The basic problem that we face is to determine the cor- 
relators for the fluctuating electromagnetic field existing 
in the empty gap between the plates. This can be done 



by suitably generalizing the methods used in heat trans- 
fer studies [3j, which are also at the basis of the recent 
out-of-equilibrium Casimir investigations [jl . Both are 
based on the well known Rytov's theory [22[ of electro- 
magnetic fluctuations. On the basis of this theory, the 
field in the gap can be interpreted as the result of multiple 
scatterings off the plates surfaces, of the radiation fields 
originating from quantum and thermal fluctuating polar- 
izations within the plates. Importantly, the local charac- 
ter of the polarization fluctuations implies that the two 
plates radiate independently from each other. Therefore, 
in a stationary configuration, the radiation from either 
plate is the same as the one that would be radiated by 
that plate, if it were in equilibrium with the environment 
(at its own temperature) , the other plate being removed. 
This physical picture permits to separate the problem of 
determining the fluctuating field in the gap in two sep- 
arate steps. In the first step, one determines the field 
radiated by either plate in isolation, a problem that can 
be solved by using the general equilibrium formalism. 
The two plates are considered together only in the sec- 
ond step, where the intracavity field is finally determined, 
by taking account of the effect of multiple scatterings on 
the radiation fields radiated by the plates, as found in 
step one. 

After these general remarks, we can now start our two- 
step computation of the fluctuating e.m. field in the gap. 
We let {x,y,z} cartesian coordinates such that the vac- 
uum gap is bounded by the planes z — and z = a re- 
spectively, with plate one (two) lying at the left (right) of 
the z = (z = a) plane. We suppose that the lateral sizes 
L x and L y of the plates are both much larger than the 
separation a: L x ,L y a Boundary effects being negli- 
gible, it is mathematically convenient to impose periodic 
boundary conditions on the fields in the (x, y) directions, 
on the opposite sides of the plates. When dealing with 
the e.m. field, it is sufficient to consider the electric field 
E(t, r) only, for the magnetic field B(t, r) can be obtained 
from E(i, r) by using Maxwell Equations. The geometry 
being planar-like, the electric field in the gap can always 
be expressed as a sum of (positive-frequency) plane- wave 
modes of the form 

where 

41(-;r)=ea x He* (±) -. (2) 

Here uj is the frequency, and is the projection of the 
wave-vector onto the (x, y) plane. Periodicity in (x, y) 
directions implies that the wave-vectors kj^ belong to 
a discrete set labelled by two integers (n x ,n y ): k x — 
2irn x /L x , k y — 2itn y /L y . The index a — s,p denotes 
polarization, where s and p correspond, respectively, to 
transverse electric and transverse magnetic polarizations. 
The superscripts (+) and (-) in Eqs. |T]) and @ re- 
fer to the direction of propagation along the z-axis, the 
(+) and (— ) signs corresponding to propagation in the 
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positive and negative z directions, respectively. More- 
over, = kj_ ± k z z, where k z = ^Juj 2 jc 2 — k"]_ (the 
square root is defined such that Re(fc 2 ) > 0, Im(fc z ) > 0), 



(±) 



We 



e$» ^xki, e$» = (c/ W )k(±) x e 
note that for lu/c > kj_, when k z is real, the modes 
^ok represent propagating waves, while for uj/c < k±, 
when k z is imaginary, they describe evanescent modes. 
It is opportune to introduce a shortened index nota- 
tion, that will prove useful in the sequel. We shall use 
a roman index i to denote collectively the i-th compo- 
nent of a vector and the position r, while a greek in- 
dex a will denote collectively the polarization a and the 
wave- vector k^. In this notation the i-th component 
of £^ ( w ; r ) shall be denoted as £^\lo). Similarly, 
a kernel A a ^ ± . a iy (w) shall be denoted as A a ^ a > (uj). 
We also set XL = / dui/(2n), 5 U>U ' = 2ttS(uj — u>'), 

Eq = V^E^.n, E a > and <W' = A <5„ x t< S ny ,n'J a , a ' , 

where A = L x L y is the area of the plates. Finally, for 
any kernel A a ^ ± . a >y ± , we define 



Having set our notations, we pass now to step one. 



A. Step one: the field radiated by a single plate in 
thermal equilibrium 

As explained above, we begin by considering each 
plate in isolation to determine its radiation, and we let 
£•(<*;) , A = 1,2 the time Fourier-transform of the field 
radiated by plate A. The total radiation field £^ cq ' A ^ (<j) 
existing, respectively, to the right of plate one and to the 
left of plate two, when either plate is in equilibrium at 
temperature Ta with the environment (the other plate 
being absent) can be expressed in the form: 



,(eq;A) , 



H=£fV)+c ov ">)+£r w H, (3) 

where £^ env (u>) describes the environment radiation, 
including vacuum fluctuations and black-body radiation, 
impinging on plate A (from the right for plate one, and 
from the left for plate two), and £^ sc ' A \uj)) is the corre- 
sponding scattered radiation. These fields have the ex- 
pansions: 



,(env;A) 



(sc;A) , 



^ A) M = E4 ±) M^ ) H 



(4) 



/^^E^M^M , (5) 



£ t c ' A] M = £ 4 ±] M sty ( W ) C v) M , (6) 



where, here and in Eqs. ([7]) and (fTO)) below, the upper 

(lower) sign is for plate one (two), and S^, is the scat- 
tering matrix of plate A, for radiation impinging on the 
right (left) surface of plate one (two). It is important to 

(A) 

note that, for fixed plates orientations, the matrix S^ a , 
depends in general on the position of some fixed 
reference point chosen on plate A. If S^, is the 

scattering matrix of plate A relative to a coordinate sys- 
tem with origin at Q<- A \ then: 



S^l = e 



k (±). x U) 5(A) ik'<=F). x W 



(7) 



The amplitudes b& (u) for the environment radiation 
in Eqs. ||SJ| and ([6]) are characterized by the following 
non- vanishing well known correlators: 

(ft^M^V)) = —^f- F(lo, Ta) Re (±\ <W<W 

(8) 

where F(uj,T) = (huj/2) coth(7ia;/(2fc B T)), with k B 
Boltzmann constant. The desired correlators for the am- 
plitudes ba (w) can now be determined by exploiting the 
following relation implied by the fluctuation-dissipation 
theorem: 

(£^' a) (oj)£^ a) *(oj')) = - F(u, T A ) <WIm G\ A) (u>) , 

(9) 

where g\ A%> (uj) is the dyadic retarded Green function of 
plate A. In the vacuum to the right (left) of plate one 

(two), the Green function g\ A \lj) can be expressed in 

(A) 

terms of the scattering matrix ', as follows: 

r [A) (, a - r m (, a + 2lTiuj2 V f (±) q( a Mt) 1 



(10) 

where GfJ (to) is the retarded Green function in free 
space: 



(-)o(-) 



J(f')c 



(11) 



with 6(z) Heaviside step-function (9(z) = 1 for z > 0, 
#(z) = for z < 0). Here, J denotes the inversion 
operator, whose action on space-indices is defined as 
J(i) = J(i,r) = (i, — r). It is useful to define the ac- 
tion of J also on polarizations, wave-vectors and prop- 
agation directions as J(a) = J(a, kj_) = (a, — kj_) and 
J((±)) = (=f). The following relations hold 



4 ±h = Jl + a a )/2 + £^ a (l - . B )/2 , 
where s a = sign(w 2 /c 2 — fcj_) and 



?(±) 



r(T) 



(±) 
J(i)J( 



a) = 



P(a)£>(=F) 



(12) 



(13) 
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where P(a) is one (zero) for s (p) polarization. The 

reciprocity relations G\P (uj) = G\f^ (lu) satisfied by the 
Green's function, as a consequence of microscopic re- 
versibility, imply via Eq. (|10j) the following important 
Onsager's relations that must hold for any scattering ma- 
trix 

q(.A) _ K ( _^P( a )+P( a ')UA) 

Upon substituting the expression for £^ eq ' A \uj) provided 
by Eqs.pp]) into the l.h.s. of Eq. ©, and after sub- 
stituting the expression of the Green function Eqs. (fTO]) 
into the r.h.s. ofEq. ©, by making use of Eqs. flgj), (fl"2|). 
(fT3|) and (fT4|) one obtains the following expression for the 

non- vanishing correlators of the amplitudes b^ A \uj): 
(b {A) {uj) V)} = S AB ^f( Wl T A ) «W 

x (eLT 5 - s {A) ^ } s {A)t + s^t,^ - sLT } 5 (A)t ) , 

(15) 

where we collected the amplitudes ba {oS) into the 
(column) vector b^ A \uj) and we set yj(p w / cw ) = 

fe n n (pw/ew) ) where n (P Q W) = Saa , (1 + 8a y 2 &nd n (ew) = 

<5aa' (1 — s a )/2 are the projectors onto the propagating 
and evanescent sectors, respectively. Eq. (|T5|) general- 
izes the well known Kirchhoff 's law (as can be found for 
example in Q) to non-planar surfaces, and it shows that 
the fluctuating field radiated by plate A is fully deter- 
mined by its scattering matrix . We remark that for 
non-planar surfaces the matrix is non-diagonal, and 
therefore the order of the factors on the r.h.s. of Eq. (fT5l) 
must be carefully respected. Now we move to step two. 



B. Step two: determination of the intracavity field 

Without loss of generality, the intra-cavity field can be 
represented as a superposition of waves of the form: 

sm = e ] M4 +) H + fo^M^M • (16) 

The intuitive physical picture of the intra-cavity field as 
resulting from repeated scattering off the two surfaces 
of the radiation field emitted by the surfaces of the two 
plates leads to the following equations for b^^uj): 

= + SW &<-> , &<-> = 6 (2) + S (2) &<+> . (17) 
Equations (fT7)) are easily solved: 

=[/ (12) 6 (D +5 (l) C/ (21) 5 (2) ^ (lg) 



where U (AB ^> = (1 - S^S^y 1 . Together with Eq. 
(fT5)l . Eqs. (HHJ) and (flUj) completely determine the intra- 
cavity field. In particular, they determine the matrix 
(ji KK ) f or the non- vanishing correlators of the intracav- 
ity field: 

(b( K \w)b^(cj')) = 6^ KK 'K (20) 

The explicit expression of C^ KK ) in terms of S^ 1 ' and 
can be easily obtained from Eqs. (TT5|) . (fT5)) and 
(jT9|l . and it is not shown for brevity. 

III. OBSERVABLES 

The above results permit to evaluate the average of 
any observables constructed out of the intracavity field. 
Typically, the observables are symmetric bilinears of the 
electric field, of the form 

6 = X>S / rf 2 rl^(t,r)a y (r,r')^(t,r'), (21) 

ij 

where Oy(r,r') = Oji(r',r). Upon defining the matrix 

ij 

xOy(ry)££'W) (22) 
the statistical average of 6 can be written as 

(°) = 2 E E TMC^'W^] . (23) 

w>0K",-K"' 

Below we shall use this formula to determine the Casimir 
force and the power of heat transfer between the two 
plates. 



A. The Casimir force out of thermal equilibrium 

As our first example, we consider the (a;, y) integral 
of the zz components of the Maxwell stress tensor Tjy, 
that provides the total Casimir force between the plates. 
After a simple computation, one finds: 

° {KK ) [Tzz] = w ( Skk ' n(PW) + 5kj(k,) n(cw) ) • 

(24) 

Evaluation of Eq. (23]) with {KK "> given by Eq. (24}, 
leads to the following representation for the unrenormal- 
ized Casimir force Fz° noq ' ) out of thermal equilibrium: 



b (~) = 5(2)^(12) 6 (1) + £,(21) 6 (2) ; 



(19) 



F (0neq) = y- _ [ F ^ iTl )J(S^,S {2) ) 
uj>0 
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+ F(lj,T 2 )J(S^,S^)} , (25) 
where J(5 (j4) , 5 (B) ) is the quantity 



+ 



+(5 (B)t E (Pw) 5 (B) +E (ew) (S (B) +i? (fl)t S (ew)\j _ (2g) 

After we add and subtract one half of the quantity 

B = F(w, T 2 ) J(5 (1) , 5 (2) ) + F(w, T x ) J(S (2) , 

from the expression inside the square brackets on the 
r.h.s. of Eq. (J5SJ), it is easily seen that Eq. PS]) can be 
recast in the form: 



Fj 0aBq >(Ti,r a ) = 



Fi 0eq) (T 1 )+F, lue,u (r 2 



(Oeq), 



AFi ncq) (Ti,T 2 ) , 



(27) 



where 



(28) 



w>0 



and 



AFi n ^(T 1; r 2 ) = ]T ^(Ffw.T!) -F(w,T a )) 



w>0 

x [J(SW,SW)- J(SW,SW)] 
Using the identity 

"1 



where 



F(u),T) = hu 



n(u,T) 



+ n{u,T) 



exp(hw/(k B T)) - 1 ' 
Eq. (|29p can be written as: 

AF jne q)(TijT2) = |^ (n([J;Tl )_ n(cJjT2)) 

x [j(s«,s (2) )- j(s<n,s<U)] . 



(29) 



(30) 



(31) 



(32) 



On the other hand, upon substituting Eq. l[26[) into the 
r.h.s. of Eq. (|28p. after a somewhat lengthy algebraic 



manipulation, it can be seen that the quantity F s 
can be further decomposed as 



Fj 0cq )(T) = A^(T) + F^(T) . 
Here, ^4^°^ (T) denotes the divergent quantity: 



A(o)(r) = 2 y^Tr c 

— ' 1,1 



OJ>0 



(Oeq) /ja 

(33) 
(34) 



As we see, this quantity depends neither on the mate- 
rial constituting the plates nor on their distance, and we 
neglect it altogether (24|. As to the second contribution 

Fz Cq> occurring on the r.h.s. of Eq. (f3"3"[) . it has the ex- 
pression 



F je q)(T) = 2Re ^£Wl TrQ 



uj>0 



k z (u^ s^s® 



(35) 



Recalling that, according to Eq. (JTJ) , the scattering 
matrices S^ A ' depend on the mutual positions of the 
plates, it is easy to verify that the above equilibrium 
force Fz has an associated free energy F(a, T) (Fz° q ^ = 
dF(a,T)/da) equal to: 



F(a, T) = 2 Im V ^11 Tla log( i _ S (l) S (3) 



(36) 



j>0 



Eqs. (f3"5| and (f3"6"|) coincide with the equilibrium ex- 
pressions, as derived within the scattering approach 
[131 . Putting everything together, after in Eq. (|2"T1) 
we remove the divergent contribution proportional to 
A(°)(Ti) + A(°)(T 2 ), we obtain the following new exact 
expression for the renormalized Casimir force between 
the plates: 



F^(T U T 2 ) = — — V 1} z — ^+AFj noq )(T 1 ,T 2 ) . 

(37) 

Some comments are now in order. We note first of all 
that the quantities F^ cq) (T) and AFj noq) (Ti, T 2 ) are both 
finite. Indeed, as we said earlier, our expression for 
Fz Cq \T) coincides with the known scattering-approach 
expression for the equilibrium Casimir force, which has 
been shown to be finite in previous studies (l9| . As 
to A^i neq) (Ti,T 2 ), it is apparent from Eq. J32]) that 
this quantity is finite, thanks to the Boltzmann factors 
n(uj,Ti). We also note that our result has the same gen- 
eral structure as the formula derived in Refs.Jl], for the 
simpler case of two plane-parallel plates. Analogously to 
that case, we indeed see from Eq. (|37|) that the non- 
equilibrium force is the sum of the average of the equi- 
librium forces, for the temperatures T\ and T 2 , plus a 
contribution AFj neq) (Ti, T 2 ), that vanishes for Tj = T 2 



p( c q) 
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(see Eq. (|32[) ) . Moreover, it is interesting to observe 
that even for T\ ^ T2, the quantity AFz ncq \ being an- 
tisymmetric in the scattering matrices of the two plates, 
vanishes if the two plates have identical scattering ma- 
trices. Such a case is realized, for example, if the two 
plates are made of the same material and if their profiles 
are specularly symmetric with respect to the (x, y) plane. 
When this happens, the non-equilibrium Casimir force is 
just the average of the equilibrium forces, for the two 
temperatures of the plates. An analogous statement can 
be found in the third of Refs. We can easily verify 
that in the case of plane-parallel homogeneous dielectric 
plates our general formula Eq. (|3T[) reproduces the result 
of Refs. [4]. In the flat case the scattering matrices of the 
plates are diagonal, and can be taken to be of the form 



_ x 

9 (2) - S ,7? (2) e 2lk * a 



(38) 



B. Power of heat transfer 

We consider now the total power W of heat transfer be- 
tween the plates. This requires that we evaluate the sta- 
tistical average of the (x, y) integral of the z-component 
S z of the Poynting vector in the gap between the plates. 
A simple computation shows that: 



0< 



KK 



\s z ] 



c 2 k z 



(-i) K (i„.n w +j« (r) n( 



(41) 

When this expression is plugged into Eq. (|23|) we obtain: 
W= ^[F(cj,T 1 )i/(5 (1) ,5 (2) )- J F 1 (w,T 2 )iJ(5 (2) ,5 (1) )] , 

(42) 

where H(S^ A \S^ B >) is the quantity 



(A) 

where R a denote the familiar Fresnel reflection coef- 
ficients. When these diagonal scattering matrices are 
plugged into Eq. ([55]). one obtains: 



F (eq) ( T ) = 4 .4 Re E E 



j>0 



i?j 2) e 2 ^° 
1 - R% ] e 2ifc *° 



-2ik r 



d(2) 



- 1 



(39) 



In the limit of infinite plates, when 

d 2 k. 



(2tt) 2 



the above formula reproduces the well known Lifshitz 
formula [lj for the Casimir force between two dielectric 
plane-parallel slabs. On the other hand, when the scat- 
tering matrices in Eq. I|38p are substituted into Eq. 
one finds: 

AFj^Ti.Ta) - A x ft£)[n(w,Ti) - n(u;,T 2 )} 

ui>0 



+ S (A) S (_ew) _ S (cw) 5 (A) t ^ jj{AB)i ^( 



(43) 



By a lengthy computation, it is possible to verify that the 
quantity H(S^\ S^ 2 ') is symmetric under the exchange 
ofS^ andS< 2 >: 



H(SW,SW)=H(S( 2 \SW) 



(44) 



By virtue of this identity, the above formula for the power 
of heat transfer can be rewritten as: 

W = »5]i 1 )[n(w 1 Ti) ~n(oj,T 2 )}H(S ( - 1 \S (2) ) . (45) 

uj>Q 

We stress once again that this formula provides an exact 
expression for W in terms of the scattering matrices of 
the surfaces. We can consider the simple special case of 
two planar slabs. When the scattering matrices for two 
planar surfaces, given in Eq. (|38p . are substituted into 
Eq. ([IS]) , the expression for the power of heat transfer 
takes the following simple form: 

W = A x h E u[n(u, Ti) - n(uj, T 2 )] 



E 



Re(fc 2 )- 



l^l 2 



lit 



(1)|2 



1-R^R^ e 2tk * a \ 2 



- 2Im(fc z )e 



— 2alm(fc z ) 



E 



2^(1 



R^\ 2 )(1-\R 



R ( a ] e 2lk * a \ 2 



Im(i?i 1) )Re(i?i 2) ) - Rc(i?i 1) )Im(i?i 2) ) 



1 



i?£ 2) e 2lk * a \ 2 



(40) 



After we substitute Eqs. (J39J) and gOJ) into Eq. (J37J), and 
upon taking the limit of infinite plates, one finds that the 
result coincides with the non-equilibrium Casimir force 
computed in Refs.[4|. 



■9(- 



-k 2 )Ae 



-2alm(fc z ) 



Im(izi x) )Im(i&' 



II 



i?£ 2) e 2ik * a \ 2 



(46) 



In the limit of large plates, the above expression coincides 
with the known formula for the power of heat transfer be- 
tween two infinite plane-parallel dielectric slabs separated 
by an empty gap 0- 
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IV. CONCLUSIONS 

In conclusion, we have developed a new exact method 
for computing Casimir forces and the power of heat 
transfer between two plates of arbitrary compositions 
and shapes at different temperatures, in vacuum. The 
method is based on a generalization to systems out of 
thermal equilibrium of the the scattering approach re- 
cently used to study the Casimir effect in non-planar ge- 
ometries [13, EH Gil • Similarly to the equilibrium case, 
we find that also out of thermal equilibrium the depen- 
dence on shape and material appears only through the 
scattering matrices of the intervening bodies. The ex- 
pressions that have been obtained are exact, and lend 
themselves to numerical or perturbative computations 
once the scattering matrices for the desired geometry are 
evaluated. Our results provide the tool for a systematic 
investigation of the shape dependence of thermal prox- 
imity effects in nanostructured surfaces, that could be 
of interest for future applications to nanotcchnology and 
to photonic crystals. In a successive publication [23j |. we 



shall use the formulae derived in this paper to compute 
the Casimir force and the power of heat transfer between 
two periodic dielectric gratings, like those considered in 
last of Refs. [l9| . The explicit form of the scattering ma- 
trices for rectangular gratings has been worked out there, 
on the basis of a suitable generalization of the Rayleigh 
expansion. At any finite order N of the Rayleigh expan- 
sion, the scattering matrices S a . a ' are of the form 

S a , a , = S(k x ,k y ) S(k x - k' x ) 6(k y - k' y ) , (47) 

where S(k x ,k y ) is a square matrix of dimension 
2(2N + 1), k x belongs to the first Brillouin zone, and k y 
is unrestricted. For scattering matrices of this form, our 
explicit formulae for the Casimir force and the power of 
heat transfer can be evaluated numerically quite easily, 
at least for sufficiently small N. 
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